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Abstract 

We introduce a new model and mathematical formulation for planning crane moves in the 
storage yard of container terminals. Our objective is to develop a tool that captures customer¬ 
centric elements, especially service time, and helps operators to manage costly relocation moves. 
Our model incorporates several practical details and provides port operators with expanded 
capabilities including planning repositioning moves in off-peak hours, controlling wait times 
of each customer as well as total service time, optimizing the number of relocations and wait 
time jointly, and optimizing simultaneously the container stacking and retrieval process. We 
also study a class of flexible service policies which allow for out-of-order retrieval. We show 
that under such flexible policies, we can decrease the number of relocations and retrieval delays 
without creating inequities. 


1 Introduction 

A container terminal is a facility where import/export/transshipment cargo containers arrive and 
are later distributed to nearby cities or loaded onto vessels for onward transportation. On the 
seaside, containers are unloaded from (or loaded to) vessels using huge berth cranes, and on the 
land side, they are stacked in (or retrieved from) a storage yard using smaller mobile cranes. Internal 
trucks transfer containers from one point to another inside the terminal, whereas external trucks 
transfer the containers from the storage yard to outside the terminal, or vice versa. 

In a typical container terminal, the storage yard is the area in which containers are stacked 
temporarily before they are loaded onto vessels or external trucks. The storage yard is usually 
divided into several bays. There are several columns in each bay, formed by stacking containers in 
tiers. The slot occupied by a container in a bay, can be specified by the tier and the column in which 
the container is stacked (see Figure [I]). In this paper, we focus on yard operations management in 
a terminal that handles import containers. 
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Upon arrival for a retrieval, the external truck driver checks in at the terminal entrance gate with 
the information of the container to be picked up, and then drives to the bay where the container 
is stacked. The container is retrieved by a yard crane and loaded onto the truck. If the container 
is not at the top of its column and is covered, the blocking containers must be relocated. These 
relocations are costly for port operators, as they are not charged to the customer. Moreover, the 
relocations slow the retrieval process, requiring external trucks to wait to receive their containers. 

In the storage area, the number of relocation moves and the average waiting time of external 
trucks are among the key performance measures (Stahlbock and Vofi [TTj). The main question, 
then, that yard managers face is as follows: given the estimated arrival and departure times of the 
containers (in a certain period of time), what is the stacking, relocating, and retrieval plan that 
will minimize the number of relocations while retrieving each container within a reasonable time? 
In this paper, we address this question by developing mathematical models and tools that capture 
several practical considerations and provide useful insights that help port operators manage yard 
crane moves. 



Figure 1: Illustration of a bay 

The above problem in a simplified setting is referred to as the Container Relocation Problem 
(CRP) or the Block Relocation Problem (BRP). The solution to the CRP is the sequence of moves 
to retrieve all containers with a minimum number of relocations, assuming that : (1) Containers 
are retrieved in a pre-defined order; (2) No container is stacked during the retrieval process; and 
(3) Only the containers that are on top of the target container are relocated. Even with this 
simplified setting, the problem is shown to be NP-hard by Caserta et al. j2], which explains why 
many proposed solution methods have been based on heuristics rather than exact optimization 
methods. Several authors consider heuristics based upon a set of rules designed to minimize the 
expected number of additional relocations as in the paper by Kim and Hong |6j, or to postpone 
(or avoid) future relocations as much as possible as presented in the paper by Unliiyurt and Aydm 
[1CZJ and in [2], Extending the objective function to the weighted sum of crane working time and 
number of relocations, Lee and Lee [7] present a three-phase heuristic that generates an initial 
feasible sequence of moves in the first phase and reduces the number of moves and crane’s working 
time in phase two and three. In other papers, authors use local search methods. Forster and 


2 









Bortfeldt 0 use a tree search and Caserta et al. [3] use the corridor method to find the best slot for 
relocating a container. Alternatively, a number of papers are based on mathematical programming 
formulations for the CRP and attempt to solve them exactly. In one of the very first papers, Kim 
and Hong [B] present a formulation and provide a branch and bound algorithm for obtaining the 
optimal locations for relocated containers. In [2j, the authors identify the optimum sequence of 
relocation moves using a binary Integer Programming (IP) formulation. Petering and Hussein [9] 
solve the same problem using a mixed IP formulation, and show that it has fewer decision variables 
and better runtime performance. Wan et al. m and Hakan Akyiiz and Lee [ 5 ] also present a 
binary IP for the CRP and then propose IP-based heuristics. For a classification scheme and a 
comprehensive literature review of the CRP and its variants in other practical applications, see the 
recent paper by Lehnfeld and Knust [Sj. 

The critical assumption behind the CRP is that the containers are retrieved in a predefined 
order, and the common aspect of all proposed mathematical models and heuristics is that they 
focus only on minimizing container relocations, and do not capture truck wait times or crane idle 
times when no truck is awaiting a retrieval. This can be significant in that such idle time allows for 
container repositioning that will result in increased timeliness of service (reduced wait times) when 
trucks do arrive. Moreover, restricting the relocation only to containers above the target container 
can lead to suboptimal solutions. 

In our work, we introduce a new mathematical model that incorporates practical details, and 
captures customer-centric service elements, such as waiting times. In our model we include service 
times (estimated departure and arrival times of containers) and determine the time of each move 
as part of the solution. Furthermore, we require that each retrieval and stacking task be completed 
within a service time window, which can be set by the port operator based on task priorities and 
acceptable wait times. In addition, our model allows for taking into account idle time and planning 
repositioning moves that can be performed during such time. We refer to this generalized model 
as the time-based CRP and present an Integer Programming (IP) formulation for it. 

One of our main contributions is a new IP formulation that jointly minimizes the number of 
relocation moves (including repositioning moves) and service time, and provides the sequence of 
moves for retrieving containers for external trucks based on various service policies. For example, 
our model can be used by port operators to minimize external truck wait time, or to give different 
waiting time guarantees to different customers to reflect relative priorities. Moreover, we expand the 
capabilities of existing approaches by determining repositioning moves that reduce overall service 
time, rather than the number of container relocations, thus moving from a focus of cost minimization 
for the port operators to one of maximization of customer service levels. 

In addition, we extend our model and formulation to the setting where containers are stacked and 
retrieved continuously, which we refer to as the Dynamic Container Relocation Problem (DCRP) 
[5J. A special case of the DCRP is to find the best slots for stacking incoming containers such 
that they can be retrieved in the future with a minimum number of relocations. The problem has 
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been addressed in a limited number of papers. Borgman et al. [lj use a heuristic to find the best 
stacking location based on the estimated departure time and the distance from the exit point. In 
another paper, Wan et al. p3] apply index-based and IP-based heuristics to the DCRP, assuming 
a fixed stacking and retrieval order. In more recent work, Hakan Akyiiz and Lee [5] propose an IP 
formulation and different heuristics that retrieve and stack containers in a pre-defined order. Also, 
some authors such as Rei and Pedroso [10] study variants of the DCRP by relaxing some of the 
constraints such as the columns height limit. 

In our model, we relax the strict stacking ordering, and jointly optimize the ordering and the 
slot for stacking an incoming container. In the stacking process, we impose pre-defined orderings 
among groups of containers (rather than for each container). The rationale for this is that typically, 
a group of containers is discharged from the vessel and is available for stacking with no difference 
in priority among the containers in the group. 

We also propose and study a new class of retrieval policies. We relax the first-come-first- 
served (FCFS) policy and allow for out-of-order retrievals. We show this flexible retrieval policy is 
surprisingly efficient and aligned with the benefits of port operators and customers. We demonstrate 
the impacts of this policy by measuring the number of relocations and service times. We prove 
and verify through experiments that such a policy results in fewer relocations and reduced service 
times. Moreover, our experiments show that flexible retrieval planning results in equitable customer 
services, an intuitively clear notion that which we will define in our context. 

The structure of this paper is as follows. In Section [2j we present our time-based model and 
formulation and present an extended formulation for the DCRP. In Section [3j we propose a flexible 
policy for retrieving containers and present the results of numerical experiments using this policy. 
Finally in Section [4] we summarize results, and suggest directions for future work. 


2 Mathematical Model and Formulation for the Time-Based CRP 


In this section, we present the mathematical model for the time-based CRP. We also develop an 
IP formulation that expands on the BRP — I formulation presented in |2]. We illustrate cases that 


arise operationally and can be handled only by our time-based model. Moreover in Section 2.1 


we give an extension of the model that is capable of handling retrieval and stacking processes and 
thus can be used for solving the dynamic container relocation problem. The mathematical model 
and formulation in this section require that containers be retrieved in a pre-defined order (FCFS 
policy). In Section [3j we relax this assumption and introduce a class of policies that allow for 
out-of-order retrievals. Table [l] summarizes the notations that will be used in the rest of the paper 
for describing the bay configuration and moves, defined as stacking, relocation, or retrieving one 
container. 

In the time-based CRP, we are given a bay with C columns and P tiers, where containers are 
stacked on top of each other up to height P (typically four or five) in each column. In a given 
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time horizon (say a day), N containers need to be retrieved from this bay and delivered to external 
trucks. The departure time of each container is the earliest time that it can be retrieved, and is 
the same as the arrival time of its corresponding truck. Throughout this paper, departure times 
are assumed to be known beforehand. 



Notation 

Description 

Bay configuration 

CxP 

A bay with C columns and P tiers. 

Container 

(C n , d n ) 

A container with label c n e {ci,C 2 , C 3 ,...}. The time-step 
at which the external truck corresponding with container c n 
arrives, is denoted by d n £ (1,2,3,...). 

Slot 

U,j] 

The slot in the r 1 column and f l tier of the bay where 

1 < i < C and 1 < j < P. Tiers are indexed from the bottom 
such that the lowest and topmost tiers are represented 
by j = 1 and j = P, respectively. 

Relocation move 

[i,j] -» [k,l] 

A relocation from slot [i, j] to slot [k, l]. 

Retrieval move 

V, j] -» out 

Retrieving a container from slot [i,j]. 

Schedule 

S € R n 

The arrival time-steps of external trucks that pick up 

containers C\, Ci ,..., Cjv. 


Table 1: Notations 

Our time-based model relies on actual service times rather than the departure orders of con¬ 
tainers. To define arrival times of external trucks for input to the model, we discretize the planning 
horizon into time-steps, where one time-step is the minimum time needed to complete one move. 
We then translate the arrival time of the trucks to time-steps. For example if the planning horizon 
is from 9 am to limn and each move takes on average 4 minutes to complete, there will be 30 
time-steps in the planning horizon. In this example if an external truck arrives at 9:15, its arrival 
time-step is 4 (note that if a truck arrives in the middle of one time-step, we assign it to the next 
time-step; so two or more external trucks may have the same arrival time, but we require that they 
be served in the order that they arrived). For describing the bay configuration and moves, we use 
the same variables as those in BRP — I [2]. However, unlike BRP — I in which the time index, t, is 
merely for determining the order of the moves, in our model it represents the actual time of the day 
that each move is performed. Moreover, the input of the model is the arrival time of the external 
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trucks rather than the departure order of containers. The four sets of variables are as follows: 


, 1 

[ 1 

if container n 

is in [i,j\ at time f, 

(1) 

Uijnt ~ j 


otherwise; 

1 

[ 1 

if container n 

is relocated from [i.j] to [k,l\ at time t, 

(2) 

ijklnt ~ 1 


otherwise; 

j 

f 1 

if container n 

is retrieved from [i.j] at time t, 

(3) 

Vijnt ~ j 

[o 

otherwise; 

Vnt = | 

[ 1 

if container n 

has been retrieved at time t €{l,...,f-l}, 

(4) 

1° 

otherwise; 



Given the arrival time of the external trucks, we can specify the time-step at which each con¬ 
tainer is ready to be retrieved. We denote the departure time of container c n (i.e., the time-step 
that c n is ready to be retrieved) by d n . Further, we denote by 5 n the maximum amount of delay 
(measured in time-step) that is allowed in retrieving container c n . Note that retrieval delay is the 
same as the wait time of external truck, and we will use these two terms interchangeably. Given 
d n and 5 n , we define a retrieval time window for each container as [d n ,d n + 5 n ], and require that 
containers be retrieved within these time windows. Also, if the departure time of a container is 
unknown or beyond the planning horizon, we allow it to stay in the bay arbitrarily long by setting 
its departure time to infinity. 

Let T be the total number of time-steps required to retrieve the containers. Because each 
retrieval task has a due time determined by its service time window, T can be expressed as follows: 


T = max {d n + 5 n }. 

1 <n<N, d n <oo 


(5) 


The set of constraints that ensures the configuration is feasible at each time-step and at most 
one move is performed at each time-step is given by (f6])-(|11|) , which are similar to the constraints 
in BRP — / [2j. We add Constraints (12)-(15) to ensure the containers are retrieved in the order 
given by their indices in the bay, and within their service time windows. 

Constraints ([b]) ensure that at each time-step, each container is either in the yard, or has been 
retrieved. Constraints ([7]) and ([8]) ensure that each slot is occupied by at most one container, and 
there are no empty slots between containers in a column. Constraints ([£]) require that at each time- 
step, at most one move takes place. Constraint (10) update the configuration of the yard at time- 
step t based on the configuration and the move that took place in time-step t — 1. Constraints 
capture the relations between retrieval variables and their corresponding configuration variables. 

Constraints (12) ensure that at most n — 1 retrievals are performed before retrieving each 
container, ( c n ,d n ). (Proposition 3.1 in Section [3] shows that satisfying Constraints (12) imply that 
containers are retrieved in prescribed order). Constraints (13) ensure that each container is retrieved 
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only after the arrival of its corresponding external truck. Constraints (14) ensure that containers 


are retrieved within their allowable time windows, and Constraints (15) require that containers are 
not retrieved at some time beyond their allowed time window. 


C P 


J^^bijnt + Vnt = 1, n = 1, t = 1,T (6) 

*'=i ;= 1 

N 

i = 1 /-/C, j = t = (7) 

n=1 

N N 

bjjnt — bij+int, i = 1,..., C, y = 1, ...,P — 1, t = 1,T (8) 

n=l n=1 

C P N C P N 

EEE x ijklnt "t EEE 3/u'nf — 1/ t = (9) 

i ; fc=l ;,/=l n=l i=l ;=1 n=l 

CP CP 

bijnt = bijnt-1 ~ x ijklnt-l + x klijnt-1 — }/ijnt-l (19) 

fc=l /=1 fc=l Z=1 

z = l,...,C, j = l,...,P, n = 1, t = 2,...,T 

c p t -1 

= n = 1 '-' N ' t = l,...,T (li) 

;=1 y=l f' =1 

CP C P N t-l 

«-">EE 1/ i/'nt T EEEE yijnf — t 1 (12) 

i=1 ;=1 i=l ;=1 n=l t'=l 

t = 1,..., T, n = 1,..., N and d\ < c ?2 < ■ ■ ■ < dw 

C P d n -1 

y"*, 7, T" 1 , yijnt = 0/ n = l,...,N (13) 

1=1 ;=1 f=l 
C df^\-5 n 

Yj yijnt = 1/ n = l,...,N (14) 

1=1 ;=1 f=rf H 
CP T 

YY Y yijnt = 0 ' n = l,...,N (15) 

i—\ j— 1 t—d n +6/7 + 1 

The objective is to jointly minimize the relocations and delay (difference between the time 
that a container is delivered to an external truck and the time that the truck arrives). We define 
the objective function as the weighted sum of these two efficiency metrics. The total number of 
relocation moves can be computed as E,fjt=i E^/=i E^=i Ef=i x ijklnt- 

T 

Also note that the earlier we retrieve container n, the larger is the summation E v nt- Thus 

t=d n 

we use the negative of this quantity as a measure of delay in retrieving containers. We also define 
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parameters w re i and w r as the weight factors for relocations and retrieval delays, respectively. The 
weight factors are set by port operator based on the port policies and sensitivity of port operations 
to the number of relocations and total delay. The optimization model for the time-based CRP is 
as follows: 


C P N T NT 

min w re i EEEE %ijklnt ^V r EE Vnt (16) 

i,k =1 j,l=l n= 1 f=l n =1 t=d„ 

si.(|6|) - (p|) and b,x,y,v> 0 

One of the main features of the time-based model is that the arrival time of external trucks can 
be specified. This is important because the trucks do not necessarily arrive uniformly over time 
(especially during off-peak hours), resulting in idle times, i.e., the time that the yard crane is not 
busy with the retrieval process. In practice, during idle times, port operators may reposition the 
containers, in order to speed up the retrieval process in the future. 

Incorporating service times into the model, we can account for the idle times and plan reposi¬ 
tioning moves based on truck arrival schedules and bay configuration. Note that the repositioning 
might not affect the total number of relocations, but it certainly affects the service delay for future 
arriving trucks, because some of the unavoidable relocations are performed beforehand (during the 
repositioning process). Thus, a reasonable way to plan repositioning is by considering its effect on 
the number of relocations and service delay (for future truck arrivals). Through jointly minimizing 
these two objectives, our time-based model can decide whether to do repositioning during the idle 
period, and which repositioning moves are most beneficial. The decision depends on the relative 
weights of the two objectives ( w re i and w r ). If w re \ is large compared to w r , the optimal decision 
will be to avoid repositioning or perform very few repositioning moves. Conversely, if w re i is low 
compared to zv r , it might be optimal to perform some additional relocations (during the reposition¬ 
ing process) in order to further decrease service delay. In Appendix A, we illustrate this trade-off 
between relocation and service delay. 

We investigated whether it is beneficial for the port (and customers) to schedule trucks in such 
a way that containers can be repositioned regularly. To examine this policy, we study the effect of 
two types of traffic (uniform and non-uniform) on service delay and on the number of relocations. 

Consider a C X P bay with N containers with a given initial configuration. Suppose that there 
are N external trucks each picking up one container. We study two cases regarding the arrival 
schedule of the external trucks: 1. Uniform traffic: N external trucks arrive at the yard uniformly 
and continuously over time, i.e., exactly one truck arrives at each time-step, and 2. Non-uniform 
traffic: N external trucks arrive at the yard in two batches (two bursts of traffic) and there is a 
period of no-arrivals between the arrival of the two batches. We generate 500 random instances 
of a 4 X 6 bay with 18 containers, and solve each instance twice: with uniform and non-uniform 
arrival times of the external trucks. In the case of non-uniform traffic, the arrival schedule is such 


that the first 9 containers arrive at times 1 to 9, and the second batch arrives at times 20 to 28; no 
truck arrives between time 9 and 20. We then compare the average delay per container (over all 
instances) between the two cases. 

Note that the number of relocations required to retrieve the containers is the same in both 
cases. To study the impact on service delay, we set the objective function as the weighted sum 
of relocations and retrieval delay with equal weights. For each instance, the moves for retrieving 
the first batch of containers (cj to eg) are the same in both cases. In 38 out of 500 instances, 
no repositioning is performed during the off-peak hours. In other instances, the service delay of 
containers in the second batch is decreased when there is non-uniform traffic because the containers 
are repositioned, i.e., some of the relocation moves can be performed during the period with no 
truck arrivals. 

Figure [2] shows the distribution of the difference of average delay per container for the two cases 
(uniform traffic delay - non-uniform traffic delay). The red line in the box is the median and the 
blue dot is the mean of the distribution for the 500 instances. The average delay (over all 500 
instances) per container is decreased by 17% when traffic is non-uniform. 



Figure 2: Effect of traffic burst on delay per container: relative decrease in delay per container for the non-uniform 
traffic case compared to the uniform traffic case (average over 500 instances for 4x6 bay) 

The impact of different scheduling policies on service delay is especially important in managing 
appointment systems and scheduling external trucks. In terminals with appointment system, truck 
drivers can make appointments to pick up or drop off containers during specific time windows and 
such a system has been shown to be very effective in reducing the number of relocation moves and 
service delay (Zhao and Goodchild jjll]). Given bay configuration and the current arrival schedule 
of external trucks, the appointment system can be used to recommend to the next truck a pick-up 
time that results in less retrieval delay and thus less waiting time. 

The second important feature of the time-based model is the service-time window. Using time 
windows allows for more flexible constraints and retrieval policies. For instance, it is possible for 
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the port operator to choose time windows such that certain trucks (for example those registered in 
the appointment system in advance) experience less delay, or to serve some trucks (retrieve some 
containers) out-of-order with specified maximum amounts of service delay for each truck. This 
policy is discussed in detail in Section [3j In addition to modelling flexibility, using time windows 
allows for more efficient computation by shrinking the search space of the IP (because we set many 
variables to zero by restricting the containers to be retrieved within a time window). 

Note that one can set the time windows arbitrarily large without affecting optimality; i.e., any 
solution to a problem with a smaller time window is also a solution to the problem with larger time 
windows. However, for computational reasons, we strive to tighten the time windows as much as 
possible. In short, we are interested in the smallest time window that ensures problem feasibility 


and optimality. Proposition 2.2 gives a bound with such properties for the retrieval time window 
of each container (<5*), for a container relocation problem with the objective to minimize the total 
number of relocations. We use the following lemma to prove Proposition |2.2| 


Lemma 2.1. Consider a CRP with the objective to minimize the total number of relocation moves. 
There exists an optimal sequence of moves, M*, such that if at any time-step, t, there is a container 
that is ready to be retrieved, that time-step is not idle. 


Proof. See Appendix B. □ 

Proposition 2 . 2 . Consider a bay with N containers, where d\ < c ?2 < ■ ■ ■ < rfjv- Suppose the 
objective is to minimize the total number of relocation moves, R* is the optimal value when b n is 
arbitrarily large, and H > R* is an upper bound given by any heuristic (easily computable). We can 
reduce 5 n to 5* n , given by without increasing the number of relocations. 


5* n = H + (n - d n y 


(17) 


Proof. In the time interval [l,d„ + <5 n ], two types of moves (relocation and retrieval) or idle time- 
steps are possible. By definition, exactly n containers must be retrieved by the end of time d n + 5 n - 
Moreover, H provides an upper bound on the number of relocation moves (H can be computed 
using existing fast heuristics). 


To count the idle time-steps in time interval [1 ,d n + 5 n \, we use Lemma 2.1 in the following way: 
Suppose container ( c n ,d n ) is retrieved at time d n > d n . The case that d n = d n is trivial because 
5 n is zero, and any bound obtained for 5 n would be valid. If d n > d n , then we know that in the 
time interval [d n ,d n ], container c n is available to be retrieved; so there is no idle time-steps in this 
interval. Also, in the time interval [1 ,d n — 1] there are at least n — 1 time-steps where a container 
is available to be retrieved, which means there are at most d n — 1 — (n — 1 ) non-idle time-steps. 
Therefore, in the time interval [1 ,d n + 5„], there are at most d n — n idle time-steps (note that 
d n = d n + 5 n because we have a solution in which c n is retrieved at d n ). By counting the maximum 
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number of relocations, retrievals, and the idle time-steps, we get an upper bound on the retrieval 
time window as follows: 

d n + b n < n + H + ma \(d n - n, 0) => 5 n < H + (n - d n ) + . (18) 

Note that <5* given by © is the largest required time window. Also, it might be tight for some 
containers, but not for all. 


□ 


2.1 The Dynamic container relocation problem 

Up to this point we have assumed that all containers are already in the bay and we need only to 
retrieve them. In practice, the retrieval and stacking processes can overlap. In other words, at the 
same time that some containers are retrieved, some other containers need to be stacked in the bay; 
we refer to this problem as the dynamic container relocation problem (DCRP). It is not clear how 
to model the DCRP using retrieval and stacking order for the containers, because a single order 
cannot be specified for both tasks. Even defining two separate orders (for retrieval and stacking) 
does not work because first, it would not be clear which of the two tasks (retrieval or stacking) 
should be prioritized. Second, the stacking process is more flexible because typically a group of 
containers needs to be stacked and there is no priority for stacking one container before others. 
We can use this flexibility to reduce the number of relocations. In this Section, we extend the 
time-based model and formulation to the DCRP and illustrate how our model can find an optimal 
sequence of stacking and retrieval moves and take advantage of flexibility in the stacking process. 


We define movement and configuration variables for the incoming containers in (19) and (20), 


and add Constraints (21) to capture the relationship between stacking variables and their corre¬ 


sponding configuration variables. 


$ijnt 


Znt 


1 if container n is stacked in [i,j\ at time f, 

0 otherwise; 

1 if container n has been stacked at time t e{l,...,f-l}, 

0 otherwise; 


(19) 

( 20 ) 


C P t-1 

z nt = YuYuYa S 'int' ' n = l,...,N, t = 1,..., T. 

i= 1 7=1 f'=l 


( 21 ) 


We also replace Constraint ©> with Constraint \22\ to ensure that at each time-step, each 
container is either in the bay, outside the bay waiting to be stacked, or has been retrieved and 
delivered to an external truck. 
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(22) 


C P 

2j 2j bijnt + Vnt = Zj;f/ M = 1/ ■■■/N, t = 1, ..., T. 

'=i y=i 


Moreover, we define a„ as the arrival time of incoming container n (the earliest time that it is 
ready to be stacked), and a n as the maximum amount of delay that is allowed in stacking container 
n, similar to d n and 5 n respectively for the arrival time of external trucks and retrieval time window. 


Note that the maximum number of required time-steps (T) is now given by 23 


T — max {.^n 

1 <n<N, 1 <m<N, d m <oo 


(23) 


We also add Constraints (24) to (26). Constraints (24) ensure that each container is stacked 
only after its scheduled arrival time. Constraints (25) ensure that containers are stacked within 
their allowable time windows, and Constraints (26) require that containers are not stacked at 
some time beyond their allowed time window. Moreover, we add stacking moves for updating 
bay configuration and replace Constraints (10) with Constraints (27). Lastly, Constraints ([9]) are 
replaced with Constraints (28) to ensure that, at each time-step, at most one move (retrieval, 
staking, or relocation) is performed. 


C P a n - 1 

^ ^ ^ Sijnt = 0, n = 1, ...,N 

i=l ;=1 f=l 
C P 

2 ^ 2 ^ H Si ’ nt = 1/ n = i '-' n 

i=l j=l t = u n 
CP T 

2], 21 H Si i nt = °' n = 1, - ,N 

i =1 j=l t=a n +a„ +1 

CP CP 

bijnt = bijiit-l ~ y 2j x ijklnt-l + ^ 2_j x klijnt-l ~~ }/ijnt-l + s ijnt-l 

k=\ 1=1 k= 1 1=1 

i = 1,...,C, j = l,...,P, n = 1, t = 2, 

C P N C P N C P N 


TjHH xi > kM + L L L ^ nt+ T jTjTj si ’ nt - 1 


t = 1.T. 


i,k=l j,l=l n=l 


i=l j =1 ii=l 


i=l j =1 n=l 


(24) 

(25) 

(26) 

(27) 

(28) 


Note that we do not need to define any order for stacking containers or sequencing stacking and 
retrieval moves; instead, we model the stacking and retrieval processes by specifying departure and 
arrival times, and time windows. Retrieval moves are performed in the predefined order (due to 
Constraints (|12|) ). It is possible to enforce stacking moves to follow a predefined ordering as well; 
however, in the stacking process usually a batch of containers (as opposed to a single container) 


12 











is discharged from a vessel and all are available to be stacked with no difference in their order. 
This is an important modelling issue that cannot be captured by specifying stacking order. In our 
proposed time-based model, we assign the same stacking time to all containers available for stacking 


at the same time and add Constraints (29) to respect the stacking order of groups of containers. 
The actual order in which the containers within a group are stacked is enforced by minimizing the 
number of relocations and total delay for the stacking and retrieval processes: 


T T 

L z " f -L z n 't, Vn,n' si. a n > = min{a,|a; > a n }. (29) 

t=i f=i 

For the DCRP, we add stacking delay and its corresponding weight to the objective function, 
the general form of which is given as follows: 


C P N T NT NT 

Wrel EEEE %ijklnt W s EE Z nt ~ W r EE V n t- (30) 

i,k =1 j,l=l n=l t=1 n=l t=a n n=l t=d n 

In Appendix C, we illustrate a case, applying the above model to find the best slots for stacking 
containers in a 3 X 4 bay such that the sum of total delay and number of relocations is minimized. 
Complete formulation of the DCRP is provided in Appendix D. 

Throughout this section, our model and formulation assumed a first-come-first-served policy 
(similar to existing models in the literature and typical practices in container terminals). In the next 
section, we relax this assumption and modify the formulation to allow for out-of-order retrievals. 
We show that this retrieval policy helps reduce the number of relocations and retrieval delay, and 
therefore benefits both port operators and customers. 


3 Flexible Retrieval Planning 


The formulation presented in Section [2] requires that containers be retrieved in a pre-defined order, 


a first-come-first-served policy, which is enforced by Constraints (12). In this section, we examine 


the FCFS in the retrieval process and propose a more general class of service policies which allows 
for flexibility in the order of retrievals. 

Suppose that there are a number of trucks in the storage yard waiting in a queue for containers 
to be retrieved and delivered to them. If the trucks are served based on a FCFS policy, the 
containers have to be retrieved in the order dictated by the arrival order of the trucks. As an 
alternative to FCFS, we introduce the following flexible retrieval policy: Suppose a container that 
needs to be relocated during the retrieval process (because it is blocking a target container) is 
ready to depart, i.e., its corresponding external truck has arrived and is waiting somewhere in the 
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queue. An alternative to relocating that container is to retrieve it out-of-order. We refer to this 
as a flexible retrieval planning policy. To avoid inequity and customer dissatisfaction, the level of 
flexibility should be defined and controlled to limit the number of out-of-order retrievals. 

We define the level of flexibility, m, as follows: for any container, ( c n ,d n ), at most m containers 
whose departure time is greater than d n , can be retrieved before c n . For example if m = 1, before 
retrieving container ( c n ,d n ), at most one container whose departure time is greater than d n may be 
retrieved. By definition, the FCFS policy is equivalent to m = 0. 

We add the following constraints to allow for out-of-order retrievals and set the level of flexibility 
m in the model (all other constraints remain the same as before.). For a bay with N containers and 
d\ < d 2 < ■ ■ ■ < dft, Constraints (31) ensure that at most n — 1 + m retrievals are performed before 


retrieving container ( d n ,c n ). The following proposition shows that if Constraints (31) are satisfied, 
each external truck experiences at most ni out-of-order retrievals. 


CP c P N t -1 

yijnt T LLLL yijnt' — t 1/ 

i= 1 j=1 i =1 j =1 n =1 f'=l 

t = 1,..., T, n = 1,... ,N and d\ < d 2 < ■ ■ ■ < d^. 


(31) 


Proposition 3.1. For a CRP with N containers, suppose that d\ < d2 < • • • < d n. If before each 
container, ( c n ,d n ), at most (n - 1) + m containers are retrieved, then the number of out-of-order 
retrievals before each container is at most m. 


Proof. See Appendix C. □ 

Example |3.2| below compares the number of relocations, retrieval delay for each container, and 
total delay, for the two cases (with m = 0 and m = 1). 


Example 3.2. Out-of-order retrieval: Consider a 3 X 4 bay with 9 containers with the initial 
configuration given in Figure [3j The arrival schedule of the external trucks that pick up the 
containers is given by S = {1,2,3,4,5,6,7,8,9}. 





(Ce, 6) 

(c 3/ 3) 

( c 4/ 4) 

(cs,5) 

(c 2 ,2) 

(Cl/1) 

(C8,8) 

(c 7 ,7) 

(0,9) 


Figure 3: Initial configuration of the bay for Example 


3.2 


We solve this instance for two cases of m = 0 and tn = 1 and compare the sequence of moves as 


shown in Figure |4a| and 4b For simplicity, the departure times of containers are not shown in the 
sequence of moves, and the containers are indicated by their indices. 
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(a) Sequence of moves for the bay of Example 


3.2 


inflexible case (m = 0, FCFS) 



I 



(b) Sequence of moves for the bay of Example 


3.2 


flexible case (m = 1) 


Figure 4 


The first two moves for both cases are relocating container c 4 to the first column and retrieving 
container C\ from slot [3,2]. Next, there are two decisions that are made differently in the case of 
m = 0 and m = 1 : 

1. After retrieving ci, container C2, which is covered by C3, is next to be retrieved if a FCFS 
policy is used. But C3 can be retrieved and loaded onto its truck (because 2 moves have been 
performed, it is time 3). In the case of m = 0, C3 is relocated to the third column. In the case of 
m = 1 , C 3 is retrieved out-of-order (before container C 2 ). 

2. After retrieving c 4 , container c 5 is next to be retrieved in a FCFS policy. But can be 
retrieved, because 5 moves (in the case of m = 1) or 6 moves (in the case of m — 0) have already 
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been performed and it is time 6 or 7. In the case of m = 0, container C(, is relocated to the third 
column. In the case of m = 1, eg is retrieved out-of-order (before container C5). 

The results of these two different decisions are that the number of relocations, total delay, and 
individual delays (denoted by W{ in Table [ 2 ]) are less in the case of m = 1. Note that the number 
of out-of-order retrievals before each container is at most 1 (although the total number of such 
retrievals is 2 in this example). In fact, in this example, the trucks that pick up containers C\, 
C3, c 4, C6, C7, cs, and C9 do not experience any out-of-order retrievals. Another important point to 
note is that even for the trucks that experience 1 out-of-order retrieval (that is, trucks that pick 
up containers C 2 and C 5 ), their wait time decreases by 1 and 2 units, respectively. 

Table [2] summarizes the results for two levels of retrieval process flexibility. Individual delays 
( W\ to w<f) denote the wait time for trucks picking up containers C\ to C9. The total number of 
relocations and total delay drop by a significant amount (67% and 32%, respectively) when we 
allow for out-of-order retrieval. Moreover, each truck has a shorter wait time when such a policy is 
in effect. 


Total 

Total 



Individual delay 




relocations 

delay 

W\ 

w 2 

W 3 ZO4 

w 5 

w 6 

W7 

w 8 

IV9 

in — 0 3 

22 

1 

2 

2 2 

3 

3 

3 

3 

3 

m = 1 1 

9 

1 

2 

0 1 

2 

0 

1 

1 

1 


Table 2: Number of relocations and delay for two levels of flexibility in Example 


3.2 


In general, we can always obtain an improved or the same quality solution by allowing out- 
of-order retrievals. As shown in proposition |3.3[ given any feasible sequence of moves, we can 
construct another sequence of moves that has the same or fewer relocations and the same or less 
retrieval delay, if any out-of-order retrieval is possible. 

Proposition 3.3. Suppose M is a feasible sequence of moves for CRP with m > 0 out-of-order 
retrievals. Given M, we can construct a new sequence of moves, M', that is feasible for CRP with 
m + 1 out-of-order retrievals and, 


(a) The total number of relocations in M' is at most the number of relocations in M; and 

(b) The delay of all containers for M' is at most the delay in M. 


Proof. The procedure for constructing M ', similar to that illustrated in Example (3.2), is to retrieve 
a container that needs to be relocated but is ready to depart. 

Suppose the total number of relocations in M is R. For all 1 < r < R, we define t(r) as the time 
of the r th relocation, and n(r) as the index of the container that is relocated at t(r). Also, we denote 
the number of out-of-order retrievals that take place before retrieving container c n , by o(c n ). 
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To construct M' , we replace each relocation 1 < r < R with retrieval, whenever the following 
two conditions are satisfied: i. d n ^ < f(r); and ii. a(c,) < m + 1, for any container (cj,dj) such 
that d{ < d n ( r ) and Cj is still in the bay. In other words, if a relocated container is ready to 
depart (condition i), we retrieve it rather than relocate it, provided that the number of out-of-order 
retrievals for each container does not exceed m + 1 (condition ii). Suppose relocation f satisfies 
condition i and ii; thus we retrieve c n ^ at t(f). For the remaining moves in the sequence (after tf), 
we form M' by replacing any moves that involve n(f) with an idle time-step. 

Now we show that M' is a feasible solution with m + 1 out-of-order retrievals. In terms of 
retrieval times, the only difference between M and M' is that some containers are retrieved earlier; 
for such containers, condition i ensures that they are retrieved only after their departure time. For 
the resulting configuration to remain feasible at each time-step, adjustments might be necessary. 
For example, suppose that there is a move in M at t > t(f) to relocate a container to slot [i,j], 
which happens to be on top of n(f). In M', we perform the same move except that we relocate the 
container to slot U,j-1 ], because container n(f) in M' is not in the bay after time t(f). Similarly, 
some of the retrieval moves in M' will be from one slot lower. Finally, condition ii ensures that the 
number of out-of-order retrievals for each container does not exceed m + 1. 

Next, we show that statement (a) and ( b ) hold for M’: If there exist at least one f that satisfies i 
and ii, then the number of relocations is decreased by at least one. Moreover, because all containers 
in M' are retrieved earlier or at the same time as in M, the delay of each container is at most its 
delay in M. 


Finally, using Lemma (2.1), some of the idle times in M' can possibly be eliminated; thus, the 


delay of some containers in M' can be reduced below the delay in M. 


□ 


In the next section, we study computationally the effect of out-of-order retrievals on the number 
of relocation moves and retrieval delay. 


3.1 Computational Experiments for Flexible Retrieval Planning 

To evaluate the out-of-order retrieval policy, we perform numerical experiments using the integer 
program introduced in Section [2] The main goal of the computational experiments is to study 
the following questions: What is the quantitative impact of an out-of-order retrieval policy on the 
number of relocation moves and on the average retrieval delay?; how does the improvement change 
as we vary bay size?; and what is the impact of an out-of-order retrieval policy on service equity? 


Effect of an out-of-order retrieval policy on the number of relocations and on delay: 

To study the impact on the number of relocations and on delay, we solve random instances for 
five different bay sizes with 4, 5, 6 , 7, and 8 columns. For each bay size, there are 4 tiers in each 
column. The initial configuration of each random instance is such that the first 3 tiers are full and 
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there is no container in the top tier; so there are initially 12, 15, 18, 21 and 24 containers in the 
bay, respectively. For each of the 5 bay sizes, we solve the random instances with different levels of 


flexibility, setting the parameter m in constraint (31) to 0, 1, and 2. Changing the bay size and the 
level of flexibility, there are 15 different test settings. We solve 1000 random instances for each test 
setting and take the average over the instances to compute the number of relocations and average 
delay (per truck). 

Figure [5] shows the results of our numerical experiments. The bay size (the number of columns 
in the bay) is shown on the horizontal axis. The percent decrease in the total number of relocations 
is shown in Figure 5a, and the percent decrease in the average delay per truck is shown in Figure 


5b In both plots, the percent decrease is computed relative to the base (m = 0, the FCFS policy) 


for two levels of flexibilitys. 




(a) 


(b) 


Figure 5: Effect of an out-of-order retrieval policy on the total number of relocations and on average 
delay 


As can be seen in Figure 5a, 1 out-of-order retrieval reduces the number of relocation moves 
by 32% in a bay of size 4x4 and by 23% in a bay of size 4x8. Note that the decrease in the 
number of relocations depends on the bay size. The percent decrease becomes smaller as the bay 
gets larger, but with a decreasing rate. To understand why, note that an out-of-order retrieval of 
container ( c n ,d n ) has two benefits: first, we avoid relocating c n , and thus reduce the total number 
of relocations by 1. Second, we avoid future relocations that might be incurred by relocating c n (if 
we relocate c n to a “bad” column that has a container with higher priority, we need to relocate c n 
again later). When the bay gets larger, it becomes more likely that we find a “good” column (a 
column that has no container with departure time smaller than d„ ) for relocating c n . Therefore, 
the beneficial effect of out-of-order retrieval diminishes as the bay gets larger. 
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The dotted line in Figure 5a illustrates that allowing for 2 out-of-order retrievals has a larger 
effect on the number of relocations. The relative decrease is 37% and 31% for a bay of size 4x4 
and 4x8, respectively. 

Figure [5b] depicts the impact of an out-of-order retrieval policy on the average delay (per truck). 
A similar trend can be seen for the relative decrease in the average delay. Allowing for 1 out-of-order 
retrieval results in a 32% decrease for a 4 X 4 bay and a 22% decrease for a 4 X 8 bay. 


Effect of an out-of-order retrieval policy on service equity: When a FCFS is replaced 
by an out-of-order retrieval policy, truck drivers might perceive service to be unfair when another 
truck is served out-of-order before them. However, a truck can receive its container out-of-order 
or can experience an out-of-order retrieval before it is served. To formalize this, we compute the 
number of out-of-order retrievals performed before serving a particular truck as follows: 


(a) 4x4 bay; m = 1 (b) 4x4 bay; m = 2 


Figure 6: Long term histogram of out-of-order retrieval moves performed before truck r is served: 
negative numbers indicate that truck r is served out-of-order; positive numbers indicate that other 
truck(s) are served out-of-order; zero indicates that neither of these cases happen 

Figure [6] shows the distribution of the number of out-of-order retrievals for truck r, for two bay 



(c) 8x4 bay; m = 1 



0.00% 10.00% 20,00% 30.00% 40.00% 50.00% 


(d) 8x4 bay; m = 2 


i 




0.00% 10.00% 20,00% 30.00% 40.00% 50.00% 
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sizes (a small bay with 4 columns and a large bay with 8 columns), and two levels of flexibility 
(m = 1 and m = 2). The integer numbers on the vertical axis show the number of out-of-order 
retrieval moves that truck r experiences in the long term and the horizontal axis shows the frequency 
of each case happening. In a flexible retrieval process, three cases can occur when truck r is served: 

• Truck r is not served out-of-order, and does not experience any out-of-order retrieval. This 
case corresponds to number zero on the vertical line. 

• One or more trucks are served out-of-order before truck r (a truck that is served out-of- 
order and also experiences some out-of-order retrievals is also classified in this group). This case 
corresponds to the positive numbers on the vertical axis. The positive numbers indicate the number 
of trucks that are served out-of-order; note that the maximum of these positive numbers depends 
on in. For example in Figures [6a] and [6cJ only +1 is on the vertical axis because m = 1. In Figures 


• Truck r is served out-of-order. For example it is the n th truck in the line, but is served as the 
in th truck where in < n. This case corresponds to negative numbers on the vertical line. The value 
of the negative numbers indicates how much earlier truck r is served. For example the case that 
truck r is the n th truck and is served as the m th truck, corresponds to in — n. Similarly, the case 
that truck r is served as the third truck is indicated by value -2. 

As shown in Figure [6j more than 40% of the time truck r does not experience out-of-order 
service for another truck and is also not served out-of-order. 

Moreover,, the frequency that truck r is served out-of-order (negative numbers) is almost the 
same as the frequency that it experiences out-of-order retrievals (positive numbers), when m — 1. 
These experiments imply that in the long term, the service that each truck receives is not adversely 
affected by of the out-of-order retrieval policy. Further, as shown in Figure [obj the average retrieval 
delay (per truck) is decreased as a result of out-of-order retrievals. 

4 Concluding Remarks 

In the storage yard of container terminals, the import containers are retrieved from the bays to 
be delivered to the external trucks (retrieval process), and, the export containers are stacked in 
the bays after being discharged from the vessels (stacking process). Service times (departure and 
arrival times of containers, or some estimate of them) are used by port operators to plan the retrieval 
and stacking process. In this paper, we introduce a time-based model and integer programming 
formulation that exploits such information. Unlike existing models, our model incorporates service 
times and thereby enables port operators to manage repositioning and plan the container movements 
when containers depart and arrive continuously. Moreover, our model captures customer service 
elements, such as service delay, and allows for adjusting levels of service. 

We also present a flexible retrieval planing policy and studied its impacts on the number of 
relocations and retrieval delay. We showed that such a policy improves operations by decreasing 


6 b and 6d there are two positive values (+1 and +2) because in 
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(or holding unchanged) the number of relocations while the waiting time of external trucks is 
shortened or remains unaffected. 

There are many papers that have proposed and studied heuristics using order-based approaches 
for the CRP and DCRP. Similar to other optimization models for the CRP, our model is also not 
suitable for real-time decision making for large instances of the problem. Thus, efforts to design 
time-based heuristics is an interesting direction for future work. Moreover, it seems that using ser¬ 
vice times in the model is more convenient when solving the CRP with incomplete information (for 
example, when we have time windows as an estimate of departure and arrival times of containers, or 
some probabilistic distribution on the service times). Therefore, studying the CRP with incomplete 
information in a time-based framework is another important direction for future research. 
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Appendix A 

Planning reshuffling moves: Trade-off between relocation and service delay 

Consider a 3 X 3 bay with 6 containers, with the initial configuration given in Figure [7} Let 
S = (1,3,4,5,6,7) be the arrival times of the trucks that pick up containers C\ to c&. Notice that 
t — 2 is an idle time-step because the first container is not blocked and can be retrieved at t = 1. For 
this simple example, there are two solutions that result in different values for number of relocation 
moves and total delay as shown in Table [3] 





(Cl/1) 

(c 2 ,3) 

(C 5 ,6) 

(C4,5) 

(c 6 ,7) 

(C3,4) 


Figure 7: Initial configuration of the bay 

In the first solution, no move is done at t = 2, whereas in the second solution, (C 5 , 6 ) is relocated 
during the second time-step, resulting in less delay for ( 03 , 4 ) later in the retrieval process. However 
in the second solution, (C 5 , 6 ) needs to be relocated again, resulting in a total of 2 relocations. 
With the objective function being the weighted sum of the relocation moves and total delay, the 
optimal solution would depend on the relative weights that determine the trade-off between the 
two metrics. For example if we assign equal weights to both terms in the objective function, both 
solutions would be optimal as they have same objective function value. If we assign a larger weight 
to relocation than delay, the first solution (with 1 relocation and 4 delay) would be optimal. 


Solution 1 Solution 2 


t=l 

retrieve (c.\, 1) : 

[1,2] - 

-> out 

retrieve (ci,l) : [1,2] 

—> out 

t=2 

idle 


relocate (C5, 6) : [3,2] 

-» [1,2] 

t =3 

retrieve (cj_, 3 ) : 

[2,2] - 

-> out 

retrieve (C2, 3 ) : [2,2] 

—> out 

t =4 

relocate (05,6) : 

[ 3 , 2 ] - 

* [2,2] 

retrieve (03,4) : [ 3 , 1 ] 

—> out 

t =5 

retrieve (C3,4) : 

[ 3 , 1 ] - 

-> out 

relocate (C5,6) : [ 1 , 2 ] 

-> [2,2] 

t=6 

retrieve (C4, 5 ) : 

[1,1] - 

-> out 

retrieve (C4, 5 ) : [ 1 , 1 ] 

—> out 

t =7 

retrieve (C5,6) : 

[2,2] - 

-> out 

retrieve (cs,6) : [2,2] 

—» out 

t=8 

retrieve (cs, 7 ) : 

[2,1] - 

-> out 

retrieve (cs, 7 ) : [ 2 , 1 ] 

—» out 

Relocations 

1 



2 


Total delay 

4 



3 



Table 3: Two sequence of moves for the bay shown in Figure |7j if w re i = w r , both solutions are 
optimal 
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Appendix B 

Proof of Lemma 12. II 


Proof. Suppose M is an optimal sequence of moves, and let M(t) be the move of time t. Also, 
suppose there exist a time step t, and a container ( Cn,d n ) such that: i. M(t ) is idle; and ii. Cn is in 
the bay at time t. and da < t. Given M, we can construct another optimal solution, M', as follows: 

We swap the two moves, M(t ) and M(t + 1). As a result, M(t + 1) becomes idle. We repeat the 
swapping for M(t + 1) and M(f + 2), if conditions (i) and (ii) are satisfied for M(t + 1). We continue 
this process until we reach a time step that does not satisfy the two conditions. In the resulting 
sequence of moves, M ', the idle time step that was initially at t is postponed to a later time step 
(that does not satisfy conditions (i) and (ii)). Next we show that M' is feasible and optimal. 

First note that an idle time step does not make any changes in the configuration; therefore 
adding or removing or shifting an idle time does not violate any of the constraints that involve 
feasibility of configuration. For the same reason, constraints [9] that ensure at most one move is 
performed at each time step, will not be violated. Secondly, note that by postponing the idle 
time step, the retrieval time of a container in M' can only be earlier than its retrieval time in 
M. Moreover, because we check condition (ii) for each swapping, it is ensured that a container 
is retrieved only after its departure time. Thus, all retrievals in M' are within the allowable time 
windows. Lastly, the containers are retrieved in the prescribed order in M', because we only swap 
an idle time step with its next move, and such a swapping does not affect the initial retrieval order 
of containers in M. 

The resulting sequence of moves, M', is also optimal; because the only difference between M 
and M' is that the order of some moves are different, which does not affect the total number of 
relocation moves. □ 
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Appendix C 

Example for DCRP 


Consider a set of 9 containers [c \,... , 09 ) that need to be stacked in a 3 X 4 bay, which is initially 
empty. The stacking schedule is given by S s = (1,1,1,1,1,2,2,2,2), which implies that the first 
five containers have to be stacked before the last four containers, but the containers in each group 


do not have to be stacked in a particular order (this is enforced by Constraint 29). The retrieval 
process starts after the stacking process (i.e. there is no overlap between the two processes in this 
example). During the retrieval process, containers C\ to c 9 will be picked up by external trucks 
whose arrival schedule is given by S r = (15,19,18,20,16,17,21,22,23). The objective is sum of 
retrieval delay, stacking delay, and total number of relocations. Solution to this problem is the 
sequence of moves for stacking and later retrieving the containers, as well as the relocation moves. 
Table 1 summarizes the solution. 


Stacking process 


t=l 

stack C5 

out —» [1,1] 

t=2 

stack Ci 

out —> [1,2] 

t =3 

stack C4 

out —» [2,1] 

t =4 

stack C2 

out —> [2,2] 

t=6 

stack C3 

out —> [ 2 , 3 ] 

t =7 

stack C9 

out —> [ 3 , 1 ] 

t=8 

stack Cs 

out —> [ 3 , 2 ] 

t =9 

stack Cy 

out —■> [ 3 , 3 ] 

t=10 

stack Cg 

out —> [ 3 , 4 ] 


t=15 
t=16 
t=17 
t=18 
t=19 
t =20 
t =21 
t =22 
t=23 


Retrieval process 


retrieve (ci, 15 ) 

[1,2] -» 

retrieve (C5, 16 ) 

[1,1] - 

retrieve (Q,, 17 ) 

[ 3 , 4 ] -» 

retrieve (C3, 18 ) 

[ 2 , 3 ] -> 

retrieve (C2, 19 ) 

[2,2] -> 

retrieve (C4, 20 ) 

[2,1] -> 

retrieve (cy, 21) 

[ 3 , 3 ] -> 

retrieve (cs, 22) 

[ 3 , 2 ] - 

retrieve (C9, 23 ) 

[ 3 , 1 ] - 


out 

out 

out 

out 

out 

out 

out 

out 

out 


Table 4: Sequence of stacking and retrieval moves 

As shown in Table [ 4 j the optimal sequence of stacking is C5, C\, C4, Ci, C3 for the containers of 
first group and C9, Cs, Cy, C(, for containers of the second group. With this stacking plan, the retrieval 
process has no relocation move and the total retrieval delay is zero. If we had to specify a fixed 
stacking order as an input to the model, the retrieval delay and relocations could be more because 
it is not obvious what is the best stacking order (that minimizes the objective). For example, if we 
set the stacking order as C\, C 2 , C 3 , C 4 , c$, eg, eg, Cy, C&, total relocations and total retrieval delay would 
be 3 and 23, respectively. Notice that as shown in this example, the time-based model for DCRP 
can also be used for determining the best slot to stack a container given its pick up time in the 
future. 
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Appendix D 

Integer programming formulation for CRP and DCRP 


c P 

H H bijnt+Vnt = 1, n = 1 t = 1 

i =1 ;'=i 

C P 

H H bijnt + V nt = Z„u n = t = 1,T 

i=lj=l 

N 

H bijnt < 1, i = 1/ —/C, j = l,...,P, t = 1,T 

n=l 


N N 

Zj — Zi bjj+int, i = 1,C, / = 1, ...,P — 1, t = 1 ,T 

n= 1 n=l 

C P N C P N 

H H H XijkM + L L L yijnt < l, t = 1 

i,fc=l ;,/=! n=l (=1 ;=1 n=l 


CPN C P N C P N 

H H H Xijklnt + E E L J/i/'rcf ZE ZE ZE Sijnt — 1/ t — 1, ..., T 

i,k =1 ;,/=! n=l i=l ;=1 ??=1 i=l ;'=1 n=l 


CP C P N t-\ 

(t - m -n)H H yijnt + H H H H yijnt' <t-i 

i =1 ;=1 j=l ;=1 n=l f'=l 

t = 1,..., T, n = 1,... ,N and di < ^2 < • ■ • < rfjv 

T T 

H z n t > H Zn't, Vn,n' si. = a„ + 1 

t=i t=i 

CP CP 

bijnt = bij n t-\ -HH Xijklnt -1 + Zj Zj Xklijnt -1 yijnt—1 

k= 1 ;=1 fc=l Z=1 

i = 1,C, j = l,...,P, n = 1, ...,N, t = 2,...,T 


CP CP 

bijnt ~ bijnt-1 ~ H H Xijklnt-1 Zj Zj Xklijnt-1 ~ yijnt—1 ® ijnt-1 

k= 1 /=1 fc=l Z=1 

z = l,...,C, j = l,...,P, n = 1, ...,N, t = 2,..., T 


c p t -1 

= E E E 

*=1 7=1 f'=l 


C P f-1 

Znf = E E E S iint >, n = 

*=i;=it'=i 
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t = 1,—,T 


for CRP 

for DCRP 

for CRP 

for DCRP 

for DCRP 

for CRP 

for DCRP 

for CRP 

for DCRP 



C P d„- 1 


EEL •///«( = o, 

i=l j=l t= 1 

n = 

C P d n +5 n 


EE E yz/nt = i, 

/=1 ;=1 t=d n 

CP T 

n = 1 

L L L y ( ;nf = 0, 

/—1 7=1 t—dyi ~$~&n +1 

n = l,...,N 

C P fl„-l 


ILL Sijnt = 0, 

n = 1, ...,N 

x=i y=i f=i 


C P a n +a n 


ELL Sijnt = 1, 
i=l 7=1 t 

CP T 

n = 1, ...,N 

EE L s ijnt = 0, 

i=l y=1 f=i7„+a„+l 

n = 1, ...,N 


for DCRP 
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Appendix E 

Proof of Proposition |3.1| 


Proof. Proof is by contradiction. Suppose there is at least one container with more than m out- 

of-order retrievals, and let n be one of such containers with m + 1 out-of-order retrievals. Since 

there are m + 1 out-of-order retrievals and we assumed that at most (n — 1) + tn containers are 

retrieved before container ( n,d n ), we know that the number of containers whose departure time is 

smaller than d n and are retrieved before container ( n,d n ), is at most (n — 1) + m — (m + 1) = n — 2. 

Therefore, at least one container (from the set {(1, d\), (2, df), ■ ■ ■, (n - l,d„_i)}) is retrieved after 

(; n,d n ). Without loss of generality, let (n - be the container that is going to be retrieved 

after ( n,d n ). Note that at least n — 1 + m retrievals (including m + 1 out-of-order) take place before 

retrieving container (n — l,d n -i), and this contradicts our assumption that before container n — 1, 

we have at most ((n — 1) — 1) + tn = n — 2 + m retrievals. 

The condition of having at most (n — 1) + m retrievals before each container n, is enforced by 

Constraint (31). These constraints ensure that if container n is being retrieved at time f, then the 

number of retrievals performed up to time t does not exceed (n — 1) + m. Note that the first term in 

the left hand side is (t — m — n) X 1 whenever container n is being retrieved at time t and therefore 

C P N t -1 

the constraint is reduced to (t - m - n) + L L L Vijnt' < f - 1, enforcing the sum of retrievals 

i= 1 j =1 n =1 f'=l 

up to time t to be at most (n — 1) + m. When container n is not being retrieved at time t, the first 

C P N t-l 

term in the left hand side is zero and the constraint is reduced to yij n r < t — 1 which 

i =1 j= 1 n=l t'=1 

is simply a loose bound (and a redundant constraint) on the number of retrieval moves up to time 
t. □ 
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